The generalized evolution of linear bias: a tool to test gravity 
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We derive an exact analytical solution for the redshift evolution of linear and scale-independent 
bias, by solving a second order differential equation based on linear perturbation theory. This bias 
evolution model is applicable to all different types of dark energy and modified gravity models. 
We propose that the combination of the current bias evolution model with data on the bias of 
extragalactic mass tracers could provide an efficient way to discriminate between "geometrical" 
dark energy models and dark energy models that adhere to general relativity. 
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PACS numbers: 98.80.-k, 98.80.Bp, 98.65.Dx, 95.35.+d, 95.36.+X 



1. INTRODUCTION 

It is well known that the large-scale clustering pattern 
of different extragalactic mass tracers (galaxies, clusters, 
etc) trace the underlying dark matter distribution in a bi- 
ased manner [l[ . Such a biasing is assumed to be sta- 
tistical in nature; with galaxies and clusters being iden- 
tified as high peaks of an underlying, initially Gaussian, 
random density field. The linear and scale-independent 
bias factor, b, is thus defined as the ratio of the mass 
tracer overdensity to that of the underlying mass over- 
density, or equivalently as the ratio of the square root of 
the mass tracer 2-point correlation function to that of the 
underlying mass correlation function. Furthermore, the 
redshift evolution of bias, b(z), is very important in order 
to relate observations with models of structure formation 
and has been shown to be a monotonically increasing 
function of redshift. 

There are two basic families of analytic bias evolution 
models. The first, called the galaxy merging bias model, 
utilizes the halo mass function and is based on the Press- 
Schechter 0] formalism, the peak-background split @ 
and the spherical collapse model 4]. Many studies have 
compared the prediction of the merging bias model with 
numerical simulations and beyond an overall good agree- 
ment, differences have been found in the details of the 
halo bias. These differences have lead to modifications 
of the original model to include the effects of ellipsoidal 
collapse @ and to either provide new fitting bias model 
parameters Q, or new forms of the bias model fitting 
function Q or even a non-Markovian extension of the 



excursion set theory Q . 

The second family of bias evolution models assumes a 
continuous mass-tracer fluctuation field, proportional to 
that of the underlying mass, and the tracers act as "test 
particles" . In this context, the hydrodynamic equations 
of motion and linear perturbation theory are used. This 
family of models can be divided into two sub-families: 

(a) The so-called galaxy conserving bias model uses 
the continuity equation and the assumption that trac- 
ers and underline mass share the same velocity field 

Then the bias evolution is given as the 
solution of a 1st order differential equation, and Tegmark 
& Peebles [TO] derived: b(z) = 1 + (b - l)/D(z), with b 
is the bias factor at the present time and D(z) the grow- 
ing mode of density perturbations. However, this bias 
model suffers from two fundamental problems: the unbi- 
ased problem ic., the fact that an unbiased set of tracers 
at the current epoch remains always unbiased in the past, 
and the low redshift problem ie., the fact that this model 
represents correctly the bias evolution only at relatively 
low redshifts z<J 0.5 [Hj]. Note that [l4j bas extended 
this model to also include an evolving mass tracer popu- 
lation in a ACDM cosmology. 

(b) A model based on the basic equation for the evolution 
of linear density perturbations, and on the assumption 
of linear and scale-independent bias, which are used to 
derive a second order differential equation for the bias, 
the approximate solution of which provides the evolution 
of bias (see [15j and [16|). The provided solution applies 
to cosmological models, within the framework of general 
relativity, with a constant in time dark energy equation 
of state parameter (ie., quintessence or phantom). 

In this article, we extend the original Basilakos & Plio- 
nis 15[ bias evolution model to provide an exact solution 
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valid for all dark energy and modified gravity cosmolo- 
gies. This implies that the current bias evolution model 
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can be used to put constraints on dark energy models 
as well as to investigate possible departures from general 
relativity. 



2. THE EVOLUTION OF THE LINEAR 
GROWTH FACTOR 

In this section, we discuss the basic equation which 
governs the behavior of the matter perturbations on sub- 
horizon scales and within the framework of any dark en- 
ergy (hereafter DE) model, including those of modified 
gravity ("geometrical dark energy"). For these cases, a 
full analytical description can be introduced by consider- 
ing an extended Poisson equation together with the Euler 
and continuity equations. Consequently, the evolution 
equation of the matter fluctuations, for models where 
the DE fluid has a vanishing anisotropic stress and the 
matter fluid is not coupled to other matter species (see 

0,0,0,&m,H|,0), is § iven b y : 
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(2.1) 



where p m is the matter density and G e g(t) — GatF(^), 
with Gn denoting Newton's gravitational constant. 

For those cosmological models which adhere to general 
relativity, [Y(t) = 1, G e s = Gn], the above equation 
reduces to the usual time evolution equation for the mass 
density contrast [2J1 , while in the case of modified gravity 
models (see [r^THl, HHHl), we have G cS ^ G N (or 
Y(t) 7^ 1). In this context, S m (t) oc D(t), where D(t) is 
the linear growing mode (usually scaled to unity at the 
present time). Changing variables from t to a, equation 
(|2.1[) becomes: 



where 
A(a) 
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(2.2) 



and B(a) 



Y(a) (2.3) 



2a 5 E 2 (a) 

with f2 m being the density parameter at the present time 
and E(a) = H(o)/Hq is the normalized Hubble function. 

Useful expressions of the growth factor have been given 
by [24j for the ACDM cosmology. Several works have 
also derived the growth factor for w(z) =const DE mod- 
els (see [2f|,[26[,[27j]), and for the braneworld cosmology 
pjj- Also Linder & Cahn [2l| derived similar expres- 
sions for "geometrical" dark energy models in which the 
Ricci scalar varies with time, as well as for models with 
a time- varying equation of state, while for the scalar ten- 
sor and f(R) models the growth factors are provided by 
Gannouji et al. [23| and Tsujikawa et al. 122] . 



etc) can be ignored (see [l3|,[il(). Then, using linear 
perturbation theory in the context of general relativity 
[Y(t) — 1, G c ff = Gn] we obtained a second order dif- 
ferential equation which describes the evolution of the 
linear bias factor, &, between the background matter and 
the mass-tracer fluctuation field: 

y5 m + 2(S m + HS m )y + 4irG eS p m 6 m y = , (3.1) 

where y = 6—1. Below, we will prove that the above 
expression is valid for any cosmological model 1 including 
those of modified gravity, with G e g = GNY{t). Since 
we also make the same assumption, as in our original 
formulation, that the tracers and the underlying mass 
distribution share the same velocity field and thus the 
same gravity field, the above equation is valid also for 
cosmological models with a modified theory of gravity. 
Using the latter we have 

j m + Vu ~ and 5 tt + Vw ~ , (3.2) 

from which we obtain 



5 tl = 



(3.3) 



Now since we assume linear biasing, we have St r = b5 m , 
and using y = b— 1, we get that d(yS m )/dt = 0. Differen- 
tiating the latter twice, we then get: yS m + 2yS m + yS m = 
0. Solving for yS m , using the fact that yS m = —y5 m and 
eq. (|2.ip we finally obtain eq. (|3.ip . 

Transforming equation (|3.1[) from t to a, we simply 
derive the evolution equation of the function y{a) [where 
y{a) = b(a) — 1] which has some similarity with the form 
of eq. (|2.2p as expected. Indeed this is 



a y 

da 2 
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da 



B(a)y = 



(3.4) 



where f(a) is the growth rate of clustering, a 
parametrization of the linear matter perturbations, given 
by: 



/(«) 



d\nS m 
dlna 



dlnD 
cflna 



(3.5) 



where fl m (a) — J7 m a~ 3 / E 2 (a) and 7 is the growth index, 
originally introduced by Wang & Steinhardt [26]. Inte- 
grating eq. (|3.5p we obtain the growth factor for any type 
of dark energy: 



D(a) 



(3.6) 



In Basilakos & Plionis [l5[, we have provided an ap- 
proximate solution of eq. (|3.ip . using f(z) ~ 1 (which is 



3. THE GENERAL EVOLUTION OF BIAS 

In Basilakos & Plionis [IH, we assumed that for the 
evolution of the linear bias, the effects of non-linear grav- 
ity and hydrodynamics (merging, feedback mechanisms 



The current theoretical approach does not treat the possibility of 
having interactions in the dark sector. Also discussions beyond 
the linear biasing regime can be found in [28j (and references 
therein). 
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valid at relatively large redshifts) , but only in the frame- 
work of general relativity, ie., Y(i) = 1, which contains a 
quintessence (or phantom) dark energy. Here our aim 
is to provide a full analytical solution for all possible 
dark energy cosmologies that have appeared in the lit- 
erature, such as a cosmological constant A (vacuum), 
time-varying w(t) cosmologies, quintessence, k— essence, 
quartessence, vector fields, phantom, modifications of 
gravity, Chaplygin gas etc. 

Inserting now y(a) = g(a)/D(a) into eq. (|3.4p and using 
simultaneously equation 12.21 and the second equality of 
equation 13. 51 we obtain: 



a g 

da 2 



A(a) 



dg_ 

da 







(3.7) 



That is, the general solution of the latter equation is 
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c\ + c 2 



da 



a 3 E(a) 



(3.8) 



where C\ and C 2 are the integration constants. Utilizing 
now a = (1 + z)-\ b = y + 1 = (g/D) + 1, b = 6(0) 
and eq. (l3.8[) . we finally obtain the functional form which 
provides the evolution of linear bias for all possible types 
of DE models, including those of modified gravity, as: 



where 



b(z) = 1 



J(z) = 



h -I 
D(z) 



C 



D(z) 



(1 + x)dx 

W) 



(3.9) 



(3.10) 



Since different halo masses result in different values of 
bo, one should expect that the constants of integration 
C\ =bo — \ and C 2 should be functions of the mass of dark 
matter halos (see [IB]), assuming that the extragalactic 
mass tracers are hosted by a dark matter halo of a given 
mass. Note that an extension of our model for the case 
of an evolving mass tracer population (ie., including the 
effects of halo merging) is provided in appendix A. 

Finally, comparing our solution of eq. (|3.10p with that 
of the usual galaxy-conserving bias evolution model, 
b(z) = 1 + (bo — 1)/D(z), it becomes evident that the 
latter misses one of the two components of the full solu- 
tion. Furthermore, our full solution does not suffer from 
the unbiased and the low redshift problems, but more im- 
portantly, the dependence of our bias evolution model on 
the different cosmologies enters through the different be- 
havior of D(a), which is affected by 7 (see equation l3.6l) . 
and of E(a) = H(a)/H Q . 

It is interesting to mention that measuring the growth 
index could provide an efficient way to discriminate be- 
tween modified gravity models and DE models which ad- 
here to general relativity. Indeed it was theoretically 
shown that for DE models inside general relativity the 
growth index 7 is well fitted by 7gr ~ 6/11 (see [2lj| . f3(il| h 
Notice, that in the case of the braneworld model of 
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FIG. 1: The bias z-evolution for different flat cosmological 
models (upper panel) and their fractional difference with re- 
spect to the ACDM model (lower panel). The models shown 
are: CPL (solid line) with w(a) — wo + Wi(l — a) and 
7 = 0.55, concordance ACDM (dashed line) and DGP (dot- 
dashed) with w(a) = — [l + f2 m (a)] _1 and 7 = 0.68. Note that 
we use fi m = 0.27, (w ,Wi) = (-0.93,-0.38) [H, b = 1.1 
and C2 = 0.45. Finally, we also plot (dotted line) the bias 
evolution for C2 = which corresponds to that of [ill ]. 



Dvali, Gabadadze & Porrati [31| (hereafter DGP) we 
have 7 w 11/16 (see also [2lj]). Indeed, it has been pro- 
posed (see [321 ]) that an efficient avenue to constrain the 7 
parameter is by determining observationally the redshift- 
dependent linear growth of perturbations. Alternatively 
other methods have been proposed in the literature, such 
as redshift space distortions in the galaxy power spec- 
trum and the growth rate of massive galaxy clusters (see 
for example [33| and references therein). It is interesting 
to mention here that the above methods also assume a 
linear and scale-independent bias. 

An alternative approach is to use the current general- 
ized bias evolution, cosmology and 7 dependent, relation 
and high quality observational bias data to test gravity. 
Of course, the observational bias data are derived for a 
particular cosmological model, but it is an easy task to 
scale them to each tested model in a consistent manner. 
Note, that such data are already available in the litera- 
ture for the case of optical QSOs (34| . If the derived value 
of 7 shows scale or time dependence or it is inconsistent 
with 7gr, ~ 6/11, then this will be a hint that the nature 
of dark energy reflects in the physics of gravity. Such an 
analysis is in progress and will be published elsewhere. 

In order to visualize the redshift and 7 dependence of 
our bias model, we compare in Fig. 1, a few flat cosmo- 
logical models in which we impose 57 m = 0.27, bo = 1.1 
and C 2 = 0.45. In particular we consider the following 
cases: 

(a) the CPL parametrization (35[ with 7 = 0.55 (solid 
line) , 

(b) the concordance ACDM (7 = 0.55, dashed line), and 

(c) the DGP with 7 = 0.68 (dot-dashed line). 

The dotted line shows the bias evolution of Tegmark 
&; Peebles [ll| model, which is also described by our bias 
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model in the limit of Ci = 0. In the lower panel of Fig.l 
we show the fractional difference of the model bias with 
respect to that of the ACDM. 



4. CONCLUSIONS 

In this work we provide a general bias evolution model, 
based on linear perturbation theory, which is valid for all 
possible non-interacting dark energy models, including 
those of modified gravity. Thus the current generaliza- 
tion of the bias evolution model can be viewed as a neces- 
sary step and an ideal tool to test the validity of general 
relativity on cosmological scales. 

It is however important to spell out clearly which are 
the basic assumptions of our model, which are common 
also to many bias models in the literature: (a) the mass 
tracers and the underline the mass share the same veloc- 
ity/gravity field, (b) the biasing is linear on the scales of 
interest (which does not preclude being scale dependent 
on small non-linear scales), and (c) that each dark matter 
halo is populated by one extragalactic mass tracer, which 
is an assumption that enters, at the present development 
of our model, only in the comparison of our model with 
observational bias data and not in the derivation of its 
functional form. 

Acknowledgments. We thank Joe Silk for useful 
comments. Manolis Plionis acknowledges funding by 
Mexican CONACyT grant 2005-49878. 



Appendix A: BIAS EVOLUTION FOR AN 
EVOLVING MASS-TRACER POPULATION 



Here we obtain the general linear bias evolution model 
assuming that the mass-tracer population evolves with 
time according to a (1 + z)" law. We now drop the as- 
sumption used in section 3, that the mass-tracer num- 
ber density is conserved in time, by allowing a contribu- 
tion from the corresponding interactions among the mass 
tracers. We obtain again the corresponding equation 
p.ip . starting from the continuity equation and intro- 
ducing an additional time-dependent term, ^(t), which 
we associate with the effects of interactions and merging 
of the mass tracers. We also make the same assumption, 
as in our original formulation, that the tracers and the 
underlying mass distribution share the same velocity field 
(or gravity field). Then: 

S m + Vu ~ and S tl +Vm + *(t) ~ , (A.l) 
from which we obtain: 



(A.2) 



Although we do not have a fundamental theory to model 
the time-dependent ^(t) function, it appears to depend 



on the tracer number density and its logarithmic deriva- 
tive as well as on the tracer overdensity: ^(t) cx ^(n, (1 + 
Str)dlnn/dt) (see eq.10 of [l4| and appendix of Basilakos 
et al. [Ta|). 

Now, in the context of linear biasing, we have 5a — b5 m 
and utilizing b = y + 1, we find that d(yS m )/dt = —ty. 
Differentiating twice the latter we then get: yd m + 2y5 m + 
yS m = — Solving for yS m , using the fact that yS m = 
—ydm — ^ and equation (|2.1[) we arrive at the following 
expression: 

y5 m +2(6 m +H6 m )y+4icG eS p m 6 m y = -2ff*-* (A.3) 



which is the corresponding equation (|3.ip for the case of 
interactions among the tracers. 

Transforming again equation (I2.I[) from t to a, we get 



da 2 



A(c 



2/(o) 



d / + B(a)y 
da 



F(a) (A.4) 



where 



F(a) = 



2#(a) + a(d^/da) 
a 2 D(a)H(a) 



(A.5) 



Now, following the same notations (y = g/D) as in sec- 
tion 3 the above differential equation becomes: 



— +A(a)— = F(a) 
da z da 



(A, 



Integrating eq. (|A.6[) it is straightforward to estimate the 
general solution of the bias factor. This is 



g(a) = Ci+C 2 



da 



da 



a 3 E{a) J a 3 E(a) 



F(a)a 6 E(a)da . 

Using the same conditions with those provided in section 
3, the bias evolution in the redshift space takes the form 



b(z) = 1 + 



D(z) 



J( z ) , Vp( z ) 
D{z) + D(z) 



where 



Vp{z) = 



m dx J0 



F(u)E(u) 
{l + uf 



du 



(A.8) 



(A.9) 



Obviously, if the interaction among the tracers is negli- 
gible (* ~ 0) then cq. (|A.8j) boils down to eq. (f3~l~0|) as it 
should. 

Now, knowledge of the exact functional form of the 
interaction term ^(z) would provide the precise redshift 
evolution of the bias. As we have a naly tically proved 
in the appendix of Basilakos et al. [f6|, a reasonable 
approach regarding the evolutionary ^(z) term is that: 
= AH (1 + z) u , where v ~ 3. Note that the Hubble 
constant has been maintained for mathematical conve- 
nience. Inserting the latter equation, a = (1 + z)^ 1 and 
d^/da = -(l+z) 2 d^ Idz into the second term of eq. (|A.4|) 
we derive that 



F(z) = A{v - 2) 



(1 + *) 



v+2 



D(z)E(z) 



(A.10) 
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where A is a positive parameter (to be determined from 
observational data see Basilakos et al. in preparation). 
Obviously, for v > 2 the derived bias evolution becomes 
stronger than in the case of no interactions, especially at 
high redshifts, which means that due to the merging pro- 
cesses the halos (of some particular mass) correspond to 
higher peaks of the underlying density field with respect 
to equal mass halos in the non-interacting case. On the 



other hand, the v < 2 case corresponds to the destruction 
of halos of a particular mass, which results into a lower- 
rate of bias evolution with respect to the non-interacting 
case. Now, for the limiting case with v = 2 we obtain 
y p = 0, implying no contribution of the interacting term 
to the bias evolution solution, as in the case with = 0, 
which can be interpreted as the case where the destruc- 
tion and creation processes are counter-balanced. 
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